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Abstract 



In this short paper we derive a formula for the spatial persistence 
probability of the Airyi and the Airy2 processes. We then determine 
numerically a persistence coefficient for the Airyi process and its de- 



' n j pendence on the threshold. 



1 Introduction 



The Airyx and Airy 2 processes are universal processes describing the fluctu- 
ation of interfaces for stochastic growth models in the Kardar-Parisi-Zhang 
(KPZ) universality class. The persistence probability is the probability that 
a process stays positive (resp. negative), or more generally, above (resp. be- 
low) a certain threshold during a time interval [0,L]. When the process 
is stationary, it might be expected that the persistence probability decays 
exponentially in L. 

The Airy processes were obtained by studying specific models in the KPZ 
universality class [3j[8j|9j|Tl|. It was only in 2010 that in an amazing experi- 
ment with turbulent nematic liquid crystals Takeuchi and Sano [13j[T5] were 
able to verify experimentally the KPZ predictions at the level of distribution 
functions and covariances (and not only at the level of the scaling exponents). 
The agreement with the theory is very good. 

In a recent paper the same authors [H] measured, among others, the spa- 
tial persistence coefficients with respect to a threshold given by the average 
of the process. In the case of the Airy 2 process, the persistence coefficients 
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have been also measured in an off-lattice Eden model [12] and verified by a 
numerical simulation of GUE Dyson's Brownian Motion [13]. 

In this short paper we determine analytic formulas for the persistence 
probability to stay below a threshold c, both for the Airyi and the Airy2 
processes. The starting point are the two works on the continuum statis- 
tics [HUD]- Then we focus on the case of the Airy! process and determine 
the associated persistence coefficient and its dependence on the threshold c. 
This is made by using the numerical approach for computing Fredholm deter- 
minants developed by Bornemann in [2]. The advantage of looking directly 
at the limit process is that we do not have uncontrolled uncertainties com- 
ing from the finite size settings of an experimental setup or of a numerical 
simulation. 
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2 Results 

In order to state the results, let us introduce some notations. We denote by 
Ai the Airy! process and by A 2 the Airy 2 process, see the review [S] for the 
definition of these processes. For a threshold c6M and a time interval [0, L] 
with L > 0, the persistence probabilities are defined by 



where A G {A 1 ,A 2 }- 

For large L, the persistence probabilities decay exponential in L with 
persistence coefficients k± given by 



As it can be seen from Figure UJ] and Figure fl~2l below, the exponential decay 
of the persistence probabilities for the Airyi process is already observed at 
relatively small values of L, for instance already at L — 1. 

The analytic result for the persistence probabilities P_ of the Airy pro- 
cesses are the following. 



P_(A, c, L) = P(A(t) < c, < t < L) 
P + (A, c, L) = P(A(t) > c, < t < L) 



(2.1) 



P±(A, c, L) ~ C±(A, c)e~ K±{Ac)L for large L. 



(2.2) 
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Proposition 2.1. For the Airyi process we have 

c, L) = det(l - ^i,l)l»CR) ( 2 - 3 ) 
where the kernel Ki t L is given by 
K ltL (x,y) = M(\x\+y + 2c) + l [x < 0] (k 1)L (x t y + 2c)-k 1)L (-x,y + 2c)) (2.4) 
with 

K 1:L (x, y) = -L= [ dz e- {x - z)2/AL e- 2L3/3 e- L{y+z) Ai(y + z + L 2 ). (2.5) 
V 47rL Jr + 

Proposition 2.2. For £/ie j4zn/ 2 process we have 

P_(A 2 , c, L) = det(l - ^ 2 , l )l 2( r) (2.6) 
where the kernel K 2 l is given by 
I<2,l(x, y) = K Ai (x + c,y + c) 

-1[ X < ] dz dfie ip ~ c)L (p{x,fi)<f)(z,fi)K AiiL (z + c,y + c) 
Jr- Jr 

(2.7) 



with 



and 



K AitL (x,y) = (e LH ^K Ai )(x,y)= [ dX e" iA Ai(A + x)Ai(A + y) (2.8) 

Jr + 

Ai(/i)Bi(x + n) - k\(x + /z)Bi(/i) 



(/>(x,fi) 



(2.9) 



Before stating the results of the numerical evaluation of (12. 3p . let us 
resume the results cited above in the following tabld^: 





A,c 


K_(A, c) 


k+(A,c) 


Experimental [14] 


A 1; 


-0.6033 


3.2(5) 


3.0(5) 


Experimental [H] 


A 2 , 


-1.7711 


0.87(6) 


1.07(8) 


Off-lattice Eden p2] 


A 2l 


-1.7711 


0.89(4) 


0.90(2) 


GUE Dyson's Brownian Motion [I~4l 


A 2 , 


-1.7711 


0.90(6) 


0.90(8) 



1 The values for the Airyi process have to be multiplied by 2 2 / 3 because the scaling 
in [2] is such that the limit process is u i-> 2 2 / 3 _4i(2 _2 / 3 u) instead of u i-> A\iu). 
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Figure 2.1: Dependence of the k_(Al,c) as a function of c. The straight 
line is the linear approximation around the average of the Airy! process 
(c = -0.6033). The slope is -4.07. 



While the different experiments and numerical simulations for the Airy2 
process provide results that are quite close to each other, no further results 
were available for the Airy! process. We first evaluated numerically the 
Fredholm determinant for the two natural thresholds, namely the average of 
the process^], c = —0.6033, and for c = 0, with the results 



A,c 


K-(A, c) 


C-{A, c) 


.Ai, -0.6033 
Ai, 


2.91 
1.10 


0.370 
0.733 



Comparing our result to the experimental one, we see that the agreement is 
fairly good. Indeed, the relative error for the Airy! process is 10%. 

We also determined the exact values of K-(Ai,c) as a function of c for 
c G [—1,0], see Figure I2TT1 We observe that this coefficient is quite sensitive to 
the threshold c (see Table [T]for the values). For instance, in the region around 
the average of the process the tangent line has a slope of —4.07, i.e., a small 
error in the centering on the threshold leads to an error in the persistence 

2 The Airyi process is a stationary process with one-point distribution given by 
P(-4i(0) < s) = F 1 (2s) where Fi is the GOE Tracy- Widom distribution function [7], 
and Fx has an average —1.20653 [T5] . 
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coefficient 4 times larger. According to the experimental results [13J and the 
analysis of specific models [B] , the approach to large time limit is at first order 
a constant times t -1 / 3 . Therefore, one has to take this finite size correction 
in account when setting the threshold (as it was made in [H]). 



3 Proofs of the analytic results 

The starting point of our analysis are two formulas on the continuum statis- 
tics for the Airyi process [10] and for the Airy 2 process [I]. Let us start with 
the Airyi process. 

Theorem 3.1 (Theorem 4 of [ID]). It holds 

P(A(t) < g(t), < t < L) = det(l - B + A L , g e- LA B ) L 2 {R) (3.1) 

where g is a function in i? 1 ([0, L]), A is the Laplacian, B (x,y) = Ai(x + y), 
and 

— (x—y) 2 /(4L) 

A L>g (x, y) = - — P m=x , b (L)= y (b(s) < g(s), < s < L) (3.2) 

V47rL 

with b a Brownian Bridge from x at time to y at time L and with diffusion 
coefficient 2. 

To get the persistence probabilities, we have to determine the explicit 
kernel for the function g(s) = c. 

Proof of Proposition \2.1[ We have to determine a formula for the Fredholm 
determinant of 1 — B + A£ j C e~ LA B . Since the Fredholm determinant is on 
all R, we can shift the variables by c and obtain the kernel 

B (x + c,y + c) - / dzA L<c (x + c,y + c)(e- LA B )(z + c,y + c). (3.3) 

Clearly, A L ^ c (x, y) = A i0 (x — c,y — c), therefore 

(E3D = Ai{x + y + 2c) - / A Lfi {x,z){e- LA B ){z + c,y + c). (3.4) 

Jr 



By the reflection principle we have 

A Lfi (x,z) = P b (o)= x ,b(L)=z(K s ) < 0,0 < s < L) 



(3.5) 



Moreover, it is known (see e.g. the review [5]) that 

e- LA B (z + c,y + c) = e -^/^+y+2c)L A ^ z + y + 2c+I *y ( 3 . 6 ) 
Putting all together we have 

(JH = Ai(x + y + 2c) - l [x<0] (K ltL (x, y + 2c) - K 1>L (-x, y + 2c)) (3.7) 
where 

K 1>L {x,y) = -L= [ dze-^ 2 / 4L e- 2L3 / 3 e- L ^Ai(y + z + L 2 ). (3.8) 
Finally, using the identity (see below) 

_J == f e - {x - z)2/iL e- 2L3/3 e- L{y+z) Ai(y + z + L 2 )dz = Ai(x + y) (3.9) 
V A-irL Jr 

we get 

K 1>L (x, y) = Ai(x + y) - K 1>L (x, y). (3.10) 

Replacing this into (13. 7p gives the desired result (12 .4p . 

Finally, let us verify ( 13. 9p . By the integral representation of the Airy 
function, 

i />e 17r / 3 oo 

Ai(6 2 + c)e 2fe3 / 3+6c = — / dw e wi l 3+bw2 - cw , (3.11) 
for any e > 0, and a Gaussian integration we get 

-^L= [ dze-^ 2 / 4L e- 2L ^ 3 e- L ^Ai(y + z + L 2 ) 

i /-e 1?r / 3 oo 

= e -L{x+y) e LV^_ / du; = + ^ ^ ^ 

2-Kl J e -in/ 3oo 

where we used again ( 13. lip . □ 



Now we consider the Airy2 process. The analogue of Theorem 13.11 for the 
Airyi process is given by 
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Theorem 3.2 (Theorem 2 of g]). It holds 

V(Mt) < 9(t), <t < L) — det(l — # Ai + A Li9 e LHAi AT Ai ) i2(R) (3.13) 

where g is a function in H l ([0,L]), H Ai = —A + x is t/ie ^4«n/ operator, 
y) = j R dXAi(x + X)Ai(y + A) zs the Airy kernel, and 



r3/ , e -(^) 2 /(4D 

(3.14) 



-x y -x3/ 3 e f _^ --p m=xm=y _ J? (6( a ) < ff ( s ) — s 2 , < s < L) 



wi/i 6 a Brownian Bridge from x at time to y — L 2 at time L and with 
diffusion coefficient 2. 

We have to determine the kernel for the special function g(s) = c. 

Proof of Proposition \2. 6 A We have to compute the Fredholm determinant of 
1 — K Ai + A LjC e~ LHAi K Ai over L 2 (M). As in the proof of Proposition I2.1[ we 
first do a shift in the variables by c and obtain the kernel 

K M (x + c,y + c) - / dzA LiC {x + c,z + c)(e LHAi K Ai ){z + c,y + c) (3.15) 
Jr 

It is easy to verify that 

A LjC (x, y) = A Lfi (x -c,y- c)e~ Lc . (3.16) 
Therefore, the kernel becomes 

dSZSD = K Ai (x+c,y+c)-e- Lc [ dz A L>0 {x, z)(e LHAl K Ai )(z+c, y+c). (3.17) 

Jr 

Thus, the desired formula follows if we can show that 

e -(x- 2 ) 2 /(4L) 

A L>0 (x,z) = e~ Lz - L /3 TP m=x , b (L)=z- L i(b(s) < -s 2 ,0< s < L) 

= l[x,*<o] / dyue ML 0(x,//)0(z,/i). 
Jr 

(3.18) 

To this end we use another representation of the kernel A^o, that can also 
be found in [I] and that follows from (13 . 1 8[) by applying the Girsanov the- 
orem and the Feynman-Kac formula. According to this characterization, 
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^l,o( x i z ) — U (L] x, z)l[ z< o] is the solution at time t = L of the boundary 
value problem 

d t u + Haiu = for x < and t e (0, L), 

u(0; x, z) = (3.19) 
u(t; z) — for a; > 0. 

The solution of this problem can be found in [TJ eq. (40)], 

u(t;x,z) = l [x<0] / dne^(f)(x,fi)(/)(z,n). (3.20) 

Note that in jl] the boundary value problem describes the action of the 
operator A^o while our formulation considers the kernel of this operator. □ 

4 Numerical approach and results 

To apply the numerical procedure of [2] we need to have an analytic kernel, 
but the kernel in Proposition 12.11 is not analytic at x — 0. This issue can 
be fixed by rewriting the Fredholm determinant as acting on L 2 (IR) into a 
Fredholm determinant acting on L 2 (R_) © L 2 (IR + ). In this way, instead of a 
scalar non-analytic kernel we get an analytic 2x2 matrix kernel. 

There are few other issues that we had to deal with by trying to compute 
numerically the Fredholm determinant with kernel (12.41) : 

• We need to introduce a cut-off T and compute the Fredholm determi- 
nant on L 2 ([—T, T]). We controlled the value of T so that by varying 
it, the result were not changing. 

• One can see that the kernel (12.41) is not bounded, but this is not a rel- 
evant problem because the conjugated kernel obtained by multiplying 
(T23D with e L{y - x) is bounded. 

• The main problem is that, even after conjugation, there are regions with 
the magnitude of the kernel that grows like e aL for a of order 1, while 
the kernel for positive x has oscillations of order 1. Consequently, the 
numerical approach works only for L relatively small because of the 
limitation due to machine precision. For the range of c G [—1,0], it 
works well at least until L = 2.5 (for c = also L = 3.5 is still fine). 

Fortunately, for the Airyi process, the logarithm of the persistence probabil- 
ity becomes rapidly a straight line as can be seen in Figure HJ] and Figure fl~2l 



S 



0.001 
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Figure 4.1: Persistence probability for the Airy! process and exponential 
interpolation ([22]) with k_ (-0.6033) = 2.91 and C_ (-0.6033) = 0.370, 
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Figure 4.2: Persistence probability for the Airyi and exponential interpola- 
tion (Q with k_(0) = 1.10 and C_(0) = 0.733. 
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below. This allowed us to determine the persistence coefficient k_ for the 
Airyx process reliably. 

Finally, let us resume in a table the values of k~(Ai, c) for c G [—1,0]. 



c 


-1.00 


-0.98 


-0.96 


-0.94 


-0.92 


-0.90 


-0.88 


-0.86 


-0.84 


-0.82 




4.858 


4.739 


4.626 


4.513 


4.402 


4.293 


4.187 


4.082 


3.978 


3.877 


C 


-0.80 


-0.78 


-0.76 


-0.74 


-0.72 


-0.70 


-0.68 


-0.66 


-0.64 


-0.62 


K- 


3.778 


3.680 


3.584 


3.490 


3.398 


3.307 


3.218 


3.131 


3.045 


2.961 


C 


-0.60 


-0.58 


-0.56 


-0.54 


-0.52 


-0.50 


-0.48 


-0.46 


-0.44 


-0.42 




2.879 


2.799 


2.720 


2.642 


2.567 


2.493 


2.420 


2.349 


2.279 


2.211 


C 


-0.40 


-0.38 


-0.36 


-0.34 


-0.32 


-0.30 


-0.28 


-0.26 


-0.24 


-0.22 




2.145 


2.080 


2.016 


1.954 


1.893 


1.834 


1.776 


1.719 


1.664 


1.610 


C 


-0.20 


-0.18 


-0.16 


-0.14 


-0.12 


-0.10 


-0.08 


-0.06 


-0.04 


-0.02 




1.558 


1.506 


1.456 


1.407 


1.360 


1.314 


1.268 


1.224 


1.181 


1.140 



Table 1: Values of K-(A.i,c) for a set of values of c G [—1,0]. The value of 
K-(Ai,0) = 1.099. 
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